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Abstract 

T production in heavy ion collisions at RHIC energy is investigated. While the transverse mo- 
mentum spectra of the ground state Y(ls) are controlled by the initial state Cronin effect, the 
excited bb states are characterized by the competition between the cold and hot nuclear matter 
effects and sensitive to the dissociation temperatures determined by the heavy quark potential. 
We emphasize that it is necessary to measure the excited heavy quark states in order to extract 
the early stage information in high energy nuclear collisions at RHIC. 
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Quarkonium production in relativistic heavy ion collisions is widely accepted as a probe of 
deconfinement phase transition at finite temperature and density. The J/iji suppression!!] has 
been observed at SPS Hi and RHICH energies and considered as a signature of the created 
new state of matter, the so-called quark-gluon plasma (QGP). At RHIC energy, however, the 
J/if/ production and suppression mechanisms are complicated, there are primordial production 
and nuclear absorption in the initial state and regeneration iQSS Si an d anomalous suppression 
during the evolution of the hot medium. 

T mesons, the bound states of bottom quarks, may offer a relatively cleaner probe into the 
hot and dense medium^ [l(| EH]. At RHIC energy, there are three important advantages in 
studying Y production compared with J /if/: (i) Y regeneration in the hot medium can be safely 
neglected[12i due to the small production cross section of bottom quarks in nucleon+nucleon 
collisions. The initial creation becomes the only production mechanism and the perturbative 
QCD calculations are more reliable for estimating the production, (ii) Ys are so heavy, there is 
almost no feed-down for them, (iii) From the Y measurement in d+Au collisions where there 
is no hot nuclear matter effect, the nuclear modification factor is about unity, implying that the 
cold nuclear absorption is negligible 113|]. Therefore, the heavy and tightly bound Ys can only 
be dissociated in the QGP phase. It is interesting to note that the recent experimental result in 
Au+Au collisions at RHIC ll411 has already indicated a strong medium effect on Y production. 

The quarkonium dissociation temperature in the hot medium can be determined by sol ving 
the Schrodinger equation for cc or bb system with potential V between the two heavy quarks 115 1 . 
The potential depends on the dissociation process in the medium. For a rapid dissociation where 
there is no heat exchange between the heavy quarks and the medium, the potential is just the 
internal energy U, while for a slow dissociation, there is enough time for the heavy quarks to 
exchange heat with the medium, and the free energy F which can be extracted from the lattice 
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calculations is taken as the potential fl6l [nil . From the thermodynamic relation F — U — TS < U 
where S is the entropy density, the surviving probability of quarkonium states with potential V = 
U is larger than that with V — F. In the literatures, a number of effective potentials in between F 



and U have been used to evaluate the charmonium evolution in QCD medium ll 1 5[ 1 1 7L 1 1 811 . 

In this Letter we investigate T mid-rapidity production at RHIC energy ( ysjvw = 200 GeV), 
by solving a classical Boltzmann equation for the phase space distribution of T states moving 
in a hydrodynamic background medium. Since the Y transverse momentum distribution should 
be more sensitive to the dynamic evolution of the system, compared with the global T yield, we 
will calculate not only the centrality dependence of the nuclear modification factor, but also its 
transverse momentum dependence and the averaged transverse momentum for the ground and 
excited T states. We will also discuss the dependence of these observables on the heavy quark 
potential and the similarity between T production at RHIC and J /if/ production at SPS. 

From the experimental data in p+p collisions, 5 1 % of the observed ground state T(l s) is from 
the direct production, and the decay contributions from the excited bb states T(lp), T(2s), T(2p) 
and T(3s) are respectively 27%, 11%, 10% and l%lfl9ll. The other two states rj b (ls) and T] b (2s) 
are scaler mesons with typical width of strong interaction and therefore not considered here. To 
simplify the numerical calculation, we will not distinguish T(2p) from T(lp) and Y(3s) from 
T(2s) and take the contribution fractions from the directly produced T(ls), T(lp) and T(2s) to 
the observed T(ls) in the final state to be ay = 51%, 37% and 12%. 

Like the description for 7/iA ll20ll . the T motion in a hot medium is characterized by the 
classical transport equation 

p"d,f T = -Cf r , (1) 

where fy(p,x,t) is the T distribution function in the phase space with T = Is, lp,2s, and the 
loss term C is responsible to T suppression in the hot medium. We have neglected here the T 
regeneration at RHIC energy, as discussed above. Taking into account the gluon dissociation 
process T + g — > b + b as the suppression source, the loss term can be written as 



where k is the gluon momentum, E g = \k\ the gluon energy, F(k, p) = k^p 11 the flux factor, 
f g (k, T, Uf/) the gluon thermal distribution as a function of the local temperature T and velocity 
11^ of the medium, and cr(k, p, T) the dissociation cross section at finite temperature. The cross 



section in vacuum can be calculated with the Operator Production Expansion method ll2 lL I22L 123 



2411 and is often used for J /if/ suppression and should be better for T suppression. The medium 
effect on the cross section is reflected in the temperature dependence of the T binding energy ey. 
By solving the Schrodinger equation for the bb system with the heavy quark potential V at finite 
temperature, one obtains ey(T) and the wave function if/(x, T) and in turn the average size of 
the system (r)(T). With increasing T, ey decreases and vanishes at the dissociation temperature 



Ty and (r) increases and goes to infinity at Ty. From the lattice simulation [25] on the J/ip 
spectral function at finite temperature, the shape of the spectral function changes only a little 
for T < T but suddenly collapses around the dissociation temperature Tj/^. To simplify our 
numerical calculation, we replace the temperature dependence of the binding energy in the cross 
section by a step function, cr(k,p,T) - <x(k, p)/&(Ty - T). Under this approximation, while 
the cross section becomes temperature independent at T < Ty, the dissociation rate a = C/Ey 
depends still on the hot medium, since the gluon density is sensitive to the temperature. 
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For V — U, H.Satz and his collaborators 1 1 5ll26ll solved the Schrodinger equation and found 
the dissociation temperatures Ty/T c = 4, 1.8, 1.6 for T = Is, Ip, 2s, respectively, where 
T c = 165 MeV|27[] is the critical temperature for the deconfinement. Since the lattice calculated 
potential is mainly in the temperature region T/T c < 4, Ty/T c - 4 is considered as a low limit 
of the dissociation temperature for the ground state. In the same way, we numerically solved 
the Schrodinger equation for V — F and the dissociation temperatures Ty/T c - 3, 1.1, 1 for 
T = Is, Ip, 2s, respectively. For all Y states, the dissociation temperatures are about 30% lower 
in case of V = F compared to that of V = U. In FigQ] we show the Y dissociation rate a as a 
function of transverse momentum at fixed temperature T c < T < Ty in the case ofV — U. In the 
calculation here we have chosen a typical medium velocity vqqp = 0.5 and assumed that vqcp 
and Y momentum p have the same direction. All the rates are large at low momentum and drop 
off at high momentum. When the temperature increases from 200 MeV (left panel) to 250 MeV 
(right panel), the rates for all Y states increase by a factor of about 2. This can be understood 
qualitatively by the temperature dependence of the gluon density, n g (T = 250 MeV)/n g (T = 200 
MeV) = (250/200) 3 = 1.95. 
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Figure 1: The T dissociation rate a as a function of transverse momentum at temperature T = 200 MeV(left panel) and 
250 MeV(right panel) for the potential V = U. The medium velocity is fixed as vqqp = 0.5 and its direction is chosen 
as the same as the T momentum. T(l.s), T(lp) and Y(2s) are respectively shown by dashed, dot-dashed and dotted lines. 
The rates for T( 1 p) and Y(2s ) are multiplied by a factor 0.1. 

In our following numerical calculations, the local temperature T{x, t) and medium velocity 
u M (x, f) which appear in the gluon distribution function f g (k, T, u) and step function ®(Ty - T) 
are controlled by the ideal hydrodynamic equations 



d^' = 0, Of, {n B if ) = 0, 



(3) 



where T^ v is the energy-momentum tensor, and the baryon density. Taking into account the 
Hubble-like expansion assumption for the longitudinal motion ll28ll . the above hydrodynamics 
describes the transverse evolution of the medium in the central rapidity region. To close the 
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hydrodynamic equations, we take the equation of state II27L 12911 for ideal parton gas with parton 



masses m u = nta — m g = and m s = 150 MeV and hadron gas with hadron masses up to 2 GeV. 
Here we did not consider the back-coupling of the Y states to the medium evolution. The initial 
condition for the hydrodynamic equations is determined by the corresponding nucleon+nucleon 
collisions and colliding nuclear geometry ll3Cill . leading to an initial temperature Tj = 340 MeV. 

The transport equation ([TJ can be solved analytically i30ll . The initial transverse momentum 
distribution fy can be described by the Monte Carlo event generator PYTHIA with MSEL=62II31 1 
for the nucleon+nucleon collisions and the Glauber model for the nuclear geometry. Since the 
collision time for the two nuclei to pass through each other at RHIC energy is less than the 
starting time of the hot medium, all the cold nuclear matter effects can be included in the initial 
distribution. Considering the fact that there is almost no Y suppression in d+Au collisions lll3ll . 
we neglect the nuclear absorption and take into account only the Cronin effect||32, 33], namely 
the gluon multi-scattering with nucleons before the two gluons fuse into an Y. The Cronin effect 
leads to a transverse momentum broadening, the averaged transverse momentum square {p])nn 
in nucleon+nucleon collisions is extended to (pj)NN + a g ^l, where / is the path length that the 
two gluons travel in the cold nuclear medium and a g ^ is a constant determined by the p+A data. 

With the known distribution function fy(p, x, f), we can calculate the yield Nj A in the final 
state by integrating fy over the hypersurface in the phase space determined by the critical tem- 
perature T c . The nuclear modification factors defined by Raa — Nj A / (NcoiiN^^ are shown in 
Figf2]as functions of the number of participant nucleons N par , in Au+Au collisions at -\fs = 200 
A GeV, where N co u is the number of binary collisions and N^ N is the Y yield in the correspond- 
ing nucleon+nucleon collisions. Since the heavy quark potential at finite temperature is not yet 
clear, we take the two limits of V = U and V = F to determine the Y dissociation temperature 
Ty. For the ground state Y(ls), the binding energy is about 1.1 GeV which is much larger than 
the typical temperature of the fireball T ~ 300 MeV at RHIC energy. Therefore, Y(ls) is un- 
likely to be destroyed in heavy ion collision at RHIC. Considering the fact that the dissociation 
temperatures Ty = 4T C for V — U and 3T C for V — F are both far above the temperature of the 
fireball, the suppression due to the gluon dissociation is very small and independent of the heavy 
quark potential. As to the excited states, the situation is different as their binding energies are 
much smaller and thus their yields are strongly suppressed. Especially, at V = F which leads to 
the lowest dissociation temperatures around the phase transition for Y(l/?) and Y(2s), almost all 
of them are eaten up by the fireball in semi-central and central collisions. Including the decay 
contribution from the excited states to the ground state, the finally observed nuclear modifica- 
tion factor for Y(ls), R AA = 2t=i.s,1p,2j fl T-VXt/ {^ coII ^'nn) * s controlled by both the ground and 
excited states. In central Au+Au collisions, as one can see in the figure, Raa ~ flij=0.51. In- 
tegrating over the number of participant nucleons, the Raa in minimum bias events is 0.53 for 
V = F and 0.63 for V = U. Both values are in reasonable agreement with the preliminary 
PHENIX data Raa < 0.6401. As one can see in the figure, while the overall Raa shows limited 
sensitivity to the heavy quark potential, the excited states display different dependence on the 
potential used in the calculation. 

In our numerical calculations, the momentum parameters (pf)j^N an d a gN are respectively 
taken to be 7.1 (GeV/c) 2 from the PYTHIA simulation JH and 0.2 (GeV/c) 2 /fm to best describe 
the available data for p+A collisions 13411 . In order to further address the Y production dynamics, 
we calculated the nuclear modification factor Raa as a function of transverse momentum p t in 
central Au+Au collisions at RHIC energy, see Figf3] There are three nuclear matter effects 
that may affect the Y p, dependence: the Cronin effect in the initial state, the leakage effect 01 
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Figure 2: Centrality dependence of the nuclear modification factors Raa in Au+Au collisions at top RHIC energy 
y^jVAi = 200 GeV. Upper and lower panels are for the heavy quark potential V = U and V = F, respectively. The 
directly produced T(l.v), T(lp), T(2.v) and the total Y(ls) are respectively shown by dashed, dot-dashed, dotted and solid 
lines. 



l30ll for higher p, particles, and the suppression mechanism in the hot medium. While both the 
Cronin effect and leakage effect lead to a p, broadening, i.e. reduction of low p, particles and 
enhancement of high p, particles, a strong enough suppression may weaken the broadening. The 
reason is the following: When the gluon traveling length / which characterizes the Cronin effect 
is long, the temperature which controls the suppression becomes high, the competition between 
the Cronin effect and the suppression may reduce the p, broadening. For the ground state, there is 
only a weak suppression, a very strong p, broadening due to the Cronin effect and leakage effect 
is expected. Since most of the excited states are dissociated in central collisions at both V — U 
and V — F, their p, broadening is strongly suppressed and not sensitive to the heavy quark 
potential. The shape of the finally observed total Raa is determined by the directly produced 
Y(ls), as shown in Figj3] 

To obtain more dynamic information on the nuclear medium through survived excited states, 
we calculated the averaged transverse momentum square (p^)aa as a function of centrality for 
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Figure 3: Transverse momentum dependence of the nuclear modification factors Raa in central Au+Au collisions at top 
RHIC energy ysjvjv = 200 GeV. The peak caused by the Cronin effect is clear at p, ~ 4.5 GeV/c for both potentials. 



Au+Au collisions, the results are shown in FigH] To reduce the theoretical uncertainty in (p})nn 
and focus on the nuclear matter effect, we considered the difference between nucleus+nucleus 
and nucleon+nucleon collisions, A(pf) = (pj)AA - (p^)nn- For the directly produced ground 
state Y(ls), the suppression is weak and the Cronin effect plays the dominant role. As a result, 
A(pf) increases monotonously with collision centrality. For the excited states, however, the effect 
of initial Cronin effect is overwhelmed by the disassociation especially in central collisions. In 
other words, in the most central collisions, the high temperature region is larger than that in 
peripheral collisions, most of the excited Ts are destroyed. Therefore, as one can see, the value 
of A(pf) goes up in peripheral collisions due to the Cronin effect, then becomes saturated in semi- 
central collisions from the competition between the Cronin effect and the increased suppression, 
and finally starts to decrease when the suppression becomes strong. In addition, as one can 
see in the figure, in case of V — F, the decrease is remarkable due to the stronger suppression 
effect. Different from the total RAA(N par t) and RAA(Pt) which are approximately one half of the 
corresponding values for the directly produced T(ls), see Figsf2]and[3] the total A(pj) is very 
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close to the one for the ground state, as shown in Fig0] 
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Figure 4: Centrality dependence of the difference A(pj) = (p^)aa ~ (p^)nn between the averaged transverse momentum 
square from Au+Au and p+p collisions at top RHIC energy yiwjv = 200 GeV. For comparison, we showed in the upper 
panel also the J/if/ data at SPS energyjj multiplied by a factor 2.4. 

It is interesting to compare the T production at RHIC energy and the J /iff production at SPS 
energy. Assuming there is no T regeneration at RHIC and no J /iff regeneration at SPS, the 
production mechanism in the two cases is then characterized by the initial creation. In addition, 
from the relationship between the fireball temperature and the dissociation temperature for the 
ground state at V = U, Ty = 4T C » Trhic and Tj^ = 2T C » Tsps, the suppression of Y(ls) at 
RHIC and J/k/j at SPS is negligible. Therefore, the averaged transverse momentum squares for 
T in Au+Au collisions at RHIC and J/^i in Pb+Pb collisions at SPS are related to each other 
through the Cronin effect, 



A(^)|f IC 



..RHIC r> 
a gN «Au 

„SPSp 

a gN «Pb 



(4) 



where we have taken a^PS = 0.08 GeV 2 /fmj32l|33|], and R Au and Rpb are the nuclear radii. This 



relation in fact predicts that the centrality dependence of A(pf) for T at RHIC is proportional to 
that for J/f at SPS. In the upper panel of Fig0]we showed the J If data at SPS|Ql multiplied 
by the factor 2.4. It is clear that the relation © works well. For RAA(N par r) and RAA(Pt), their 
behavior depends on the details of the hot medium, it becomes difficult to obtain similar relations 
between T at RHIC and J/f at SPS. 
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Figure 5: Initial temperature dependence of the nuclear modification factor Raa in central Au+Au collisions at top RHIC 
energy tJsnn = 200 GeV. 



How hot is the fireball formed in relativistic heavy ion collisions? This is a crucial question 
that will have an influence on all signatures for QGP formation. In all above calculations, the 
initial temperature T, = 340 MeV is determined by the initial energy density and baryon density 
which are controlled by the nucleon+nucleon collisions and the nuclear geometry. To extract 
the initial temperature of the system from T production, we now take Tj as a free parameter 
and calculate the momentum integrated Raa as a function of Tj in central Au+Au collisions at 
RHIC energy, the result is shown in Fig|5] While in case of V = F, see bottom plot of Figj5] 
the values of Raa are almost constants in the temperature region 1.5 < Tj/T c < 2 which is the 
expected initial temperature region for collisions at RHIC, for V = U, see top plot of FigQ the 
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excited states and the finally observed ground state are sensitive to the temperature. Therefore, 
the experimental results of Raa for any state will allow us to extract the information on the initial 
temperature of the system. 

In summary, we studied Y production in high energy nuclear collisions at RHIC in a transport 
model. The observed Y(ls) is mainly from the direct production, and the contribution from the 
feed down of the excited states is small. The transverse momentum distribution of Y(ls) is 
not sensitive to the hot medium, but characterized by the Cronin effect in the initial stage. The 
above conclusion is almost independent of the heavy quark potential. However, the behavior of 
the excited T states is controlled by the competition between the cold and hot nuclear matter 
effects and sensitive to the heavy quark potential. Therefore, the yield and transverse momentum 
distribution for the excited states should be measured in the future experiments, in order to probe 
the dynamic properties of the formed fireball. The initial state Cronin effect can be studied via the 
centrality dependence of (pf) in A+A collisions. We did not address the influence of the parton 



distribution in A+A collisions 135L 13611 . Although it will affect the details of the predictions made 
in this letter, the qualitative trends, especially the nature of the high sensitivity of the excited Y 
states to the potential and initial temperature, will remain to be true. The consideration on the 
velocity dependence of the dissociation temperature and the influence of reduced binding energy 
in low temperature region will be discussed in the future. 
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